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The theory of quantum dynamical semigroups within the mathematically rigorous framework 
of completely positive dynamical maps is reviewed. First, the axiomatic approach which deals 
with phenomenological constructions and general mathematical structures is discussed. Then ba- 
sic derivation schemes of the constructive approach including singular coupling, weak coupling and 
low density limits are presented in their higly simplified versions. Two-level system coupled to a 
heat bath, damped harmonic oscillator, models of decoherence, quantum Brownian particle and 
Bloch-Boltzmann equations are used as illustrations of the general theory. Physical and mathemat- 
ical limitations of the quantum open system theory, the validity of Markovian approximation and 
alternative approaches are discussed also. 



I. INTRODUCTION 



Classical theory of systems interacting with environment employs evolution equations which can be written in the 
following abstract mathematical form 

>. 

§§: Tt p t = c Pt , t >o. a) 



l/") | Here pt is a time dependent probability distribution over the relevant system's configuration space f2. In particular 
f2 can be either a phase-space of a system , its position or momentum (velocity) space or a discrete set in the case of 
coarse-grained description or a discretized model. The most general "continuous" form of £ is a sum of a differential 
operator (Fokker-Planck or diffusion type) and a collisional integral one [1]. For a discrete version we can write (1) 
as a Pauli master equation 

I ' 

370* (&) = ^2(akiPt(l) - ai k pt(k)) . (2) 

3 : 

Q H . The "gain — loss" structure of the RHS of (2) guarantees preservation of ^2 k Pt{k) while positivity of the transition 
probabilities per unit time (au > 0) implies P t(k) > 0. Similar mathematical consistency conditions are satisfied for 
the continuous cases due to the positivity of the diffusion constant (or positively defined diffusion matrix) and/or 
positivity of the scattering cross-sections. The solution of eq.(l) is given in terms of a Markov semigroup Tt = exp{t£} 
satisfying the composition law 



T t T s = T t+S , t, s > (3) 

and preserving positivity and normalization of the probability distribution p t . 

In the end of 60-ties and begining of 70-ties one needed an analogical formalism to describe a variety of irreversible 
phenomena in quantum optics, solid state physics, magnetic resonance, nuclear and particle physics, etc. [2,3]. The 
classical probability distribution p t should be replaced by the density matrix p t satisfying the analog of eq.(l) 

j t Pt = Lp t , t > . (4) 

Besides many heuristic attempts the mathematically sound theory of quantum dynamical semigoups has been devel- 
oped very soon [4-7]. It is convenient for the further discussion to distinguish two approaches in this theory - the 
axiomatic and the constructive one. 

The former approach, presented in the next Section, is concentrated on general mathematical properties of eq.(4) 
and its solutions which are by no means trivial and lead to the theory of comptelely positive maps and semigroups 
on operator algebras. The fundamental structural theorems will be given and illustrated by several examples. 
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In constructive approach one tries to derive eq.(4) from first principles. Starting with a model of a quantum system 
interacting with a quantum environment one obtains Markovian master equation (4) as an approximation to the exact 
reduced dynamics of the open system. The canonical models of reservoirs and basic approximation procedures are 
discussed in Section III. 

In the recent years quantum theory of open systems and in particular quantum dynamical semigroups became 
again a subject of intensive investigations. This is motivated by the new experimental techniques which allow precise 
monitoring of decohcrcnce and dissipation in microscopic and mesoscopic quantum systems [8]. Possible future 
applications of controlled quantum systems in quantum information processing and quantum computing are another 
sources of this revival [9]. Both mentioned above approaches are necessary to create physically correct, mathematically 
consistent and convenient for practical applications theory. 

The number of publications on this topic is enormous and therefore the presented list of references is far from being 
complete, often review papers, books and other contributions to this volume are quoted instead of the original papers. 



II. AXIOMATIC APPROACH 



Consider an open quantum system with the associated Hilbert space H with the time evolution governed by the 
master equation (4) for its density matrix p t . This is the Schrddinger picture of a time evolution which will be used 
in this paper. We leave as a simple exercise for the reader to translate all the presented results to the Heisenberg 
picture. 

In order to give a precise mathematical meaning to eq.(4) we first extent it to the whole Banach space T(TL) of 
trace-class operators acting on TL equipped with the norm \\o~ ||i = trfW) 1 / 2 . Then L is a linear, perhaps unbounded, 
operator defined on its domain D(L) C T{TL). The formal solution of eq.(4) is given by the the one-parameter 
semigroup A t , t > 

A ^ - E (5) 

n=0 

for a bounded operator L while for an unbouded one 

A t a = lim (1 - tL)~ n a (6) 

makes sense. The semigroup A t should satisfy the following conditions 

A t A s = A t+S - semigroup property (7) 



lim^Ai(7 = (7 — continuity (8) 
for a > , Ajcr > — positivity (9) 



for any a e T(H) , Tr (A t a) = Tr a — trace preservation . (10) 

We shall see that the particular feature of composed quantum systems, namely the existence of entangled states, 
makes us to replace the positivity condition (9) by a stronger one called complete positivity. 



A. Complete positivity 

Consider two well-separated open quantum systems with dynamics given by two dynamical maps A^ and A^ 2 ) 
respectively. Assume also that the initial state of a joined system is determined by the entangled density matrix 

P (12) ^£ftP?W 2) - (") 

3 

One can easily find examples (e.g. for 2 qubits) of positive trace preserving AW and A( 2 ) and p< 12 ) satisfying (11) 
such that A' 1 ) <g> A' 2 )// 12 ) is not positively defined [10]. The minimal condition which must be imposed on dynamical 



2 



maps to allow meaningful tensor products is complete positivity which means that for any n = 1, 2, 3, A ® I n is 
positive, where /„ is an identity map acting on n x n matrices (i.e. trivial dynamical map on n- level quantum system). 
Completely positive maps on operator algebras were studied already in the 50-ties and the celebrated Stinespring 
representation [11] leads to a general form of completely positive dynamical map called often Kraus decomposition 
[12] 

Ap = Y,W aP W* (12) 

a 

where W a are bounded operators on H satisfying ^2 a W*W a = 1. The decomposition (12) is highly nonunique, in 
particular the sum over a can be replaced by an integral. If H is n-dimensional then one can always find Kraus 
decomposition in terms of at most n 2 terms. 



B. Completely positive dynamical semigroups 

From now on by a quantum dynamical semigroup (QDS) we mean a family of maps A t ,t > satisfying conditions 
(7)-(10) with (9) strengthened to complete positivity. The celebrated result of Lindblad, Gorini, Kossakowski and 
Sudarshan [4,5] provides us with the most general form of a quantum Markovian master equation (MME) with a 
bounded generator 



/, : : -i[H, Pt ] + J2 VjPtV* V*Vj,P t } . (13) 



d 

3 3 



or in a slightly different but often used form 

j fPt = -i[H,p t ] + ^([ Vj ,p t V*} + [V jPu V*]) . (14) 

3 

The choice of bounded operators H = H* and Vj is again not unique and the sum over {j} can be replaced by an 
integral. To simplify the notation we put always h = 1 and fcs = 1 to have the same units for energy, frequency and 
temperature. 

The non-Hamiltonian part of the generator (13) displays the "gain - loss" structure similar to that of Pauli master 
equation (2) and determined by a completely positive quantum transition map 

p~$p = J2v j pv; (15) 

3 

which is a sum of "pure" transition maps p i— > VjpV* corresponding to elementary irreversible processes [7]. 

The following expansion for A t = cxp(tL) involving only sums (or integrals) and compositions of manifestly com- 
pletely positive maps $ (15) and W t 

W t p = S t pS; , S t = exp{-ztff - (t/2) V 3* V 3} ( 16 ) 

is often useful 

At = W t + Y] dt n dt„_i.../ dt{W t -t n 3>W tn -t n _^ ■ ■ .W t2 -t^W tl . (17) 
n=1 Jo Jo Jo 

For open systems with infinite dimensional Hilbert spaces we expect that the generators of QDS are typically un- 
bounded. Although there exists no theorem giving the most general structure of L in this case, the standard form 
(13) (14) makes sense very often [13,14]. Then the operators H and Vj can be unbounded, in fact Vj can be even 
operator-valued distributions (e.g. quantum fields) and the sum over {j} can be replaced by an integral. The typical 
mathematical condition which could be imposed is the existence of the contacting semigroup S t (16) on the Hilbert 
space such that the operators S t Vj became bounded and the expansion (17) makes sense. One of the unexpected 
features of the solutions to MME with unbounded generators is that despite the formal trace preservation of (13) (14) 
we can have Trp t < Trpo- This phenomenon is known for classical Pauli master equations with infinite number of 
states where for certain transition probabilities "the system can escape to infinity in finite time" [15]. 
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C. Examples 



We present few examples of QDS generators obtained using axiomatic construction based on phenomenological 
arguments. For other interesting models of dissipative quantum systems studied using numerical computations also, 
sec [16]. 

1. 2-level system 

The simplest and the oldest example of QDS for a 2-level system (qubit) can be constructed using three elementary 
transition maps. Assume that |1 > and |2 > form a basis of Hamiltonian eigenvectors and define the following 
operators 

P k = \kxk\ , fc = l,2, a+ = (O* = |2 >< 1|, cr 3 = P 2 -P 1 ■ (18) 
The transition maps correspond to three different irreversible processes 

damping p ji<t~ pa + (19) 

pumping p i— > 7jcr + ( 0(7~ (20) 

pure decoherence p 8\P\pP\ + PipPi . (21) 
The master equation obtained from (18)-(21) reads 

J t Pt = -y[^3,Pt] + ^([<T-,PtO + ] + [<7-p t ,<7+]) + ^{[cT + ,Pt<T-} + [<T + p U CJ-]) - -[<T 3 , [<T 3 , Pt]] (22) 

where 8 = (8\ +82)/^- Using the decomposition 

p t = Pl (t)P 1 + (1 - pi{t))P 2 + a(t)<j+ + a(t)a- (23) 

we obtain 

Pi(t) = Pi(0)e-^+^ t + — H_(l- e -(TT+7i)t) } «(*) = a (0) exp{-twt - (7, + 7 T + 5)i/2}. (24) 

71 + 7T 

The diagonal elements of p t evolve independently of the off-diagonal ones and the former satisfy Pauli master equation 
with a stationary state which can be written in a form of the Gibbs state 

Peq = (2cosh^)- 1 e — ^ , H = e -»/T (25) 
zl 71 

Moreover, for any initial state p t — > p eq for t — > 00 and the generator satisfies quantum detailed balance condition. 

The equation (22) written explicitly in terms of matrix elements is called Bloch equation and is frequently used in 
quantum optics and magnetic resonance theory. 

2. Harmonic oscillator 

Another well-known example is a linearly damped and pumped harmonic oscillator with a Hamiltonian H = ua + a 
([a,a + ] — 1) and transition maps describing coherent damping and pumping 

damping pi-^Jiapa + , pumping p^^/^a + pa. (26) 

The master equation with an unbounded generator 

j t Pt = -iu[a + a, Pt ] + 2L([a,p t a + ] + [ap u a + ]) + ^-([a + ,p t a] + [a + p u a]) (27) 
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is exactly solvable and the solution can be presented for instance in the form of the quantum generating function of 
the complex variables z,z [7] 

F t (z, z) = Ti(p t cxp{za - za+}) = e - A ^F (z u z t ) (28) 

where 

z t = ze X p{-iujt-( 7l -^)t/2} , A(t) = ^— !L(i_ e -(7i-7T)*) . (29 ) 

By differentiation of F t (z, z) with respect to z or z one obtains time dependence of all moments Tr(p t (a + ) n a m ). One 
can see from (28) (29) that if 7^ > 7-f then Hindoo p t — p eq where 

p eq = (1 - e-/ T ) expj-ica+a} , £ = log ^ (30) 

j j 7l 

One can also show that the diagonal and off-diagonal matrix elements of pt evolve independently with the former 
satisfying detailed balance Pauli master equation. 

Equation (27) found numerous applications, for instance in quantum optics to describe single mode of radiation 
in cavity. This is also the simplest example of the important class of the quasi-free QDS (see [7,10] and references 
therein). 



3. Decoherence of mass center 

Following [17] we briefly present a model based on QDS which can be used to analyse environmental decoherence 
which cause the transition from quantum to classical world. Consider a mass center of a body described by the 
Hamiltonian 

H = ^-P 2 + V{X) (31) 

where X and P arc the operators of mass center and total momentum satisfying 

[X k ,P l ]=i6 kl , k,l = 1,2,3. (32) 

We assume that the interaction with an environment can be reduced to processes of scattering, absorption or emission. 
A collision with a gas particle and emission, absorption or scattering of a photon (or other quasiparticle) causes a 
momentum transfer k which changes the total momentum. The following transformation (in the Heisenbcrg picture) 

e ikx pe -ikx = p + k (33) 

describes this momentum transfer independently of the detailed microscopic mechanism of energy redistribution. 
Therefore the elementary irreversible process is given by the transition map 

p ^ e- 4kx pe lkx . (34) 

Assuming statistical independence of different momentum transfer events (called simply collisions) we obtain the 
following form of the MME 

j t p t = -i[H,p t ] + J d 3 kn(k)(V* kx p te * kx - Pt ) (35) 

where n(k) is a density of collisions per unit time leading to the momentum transfer k. 

The generator given by (35) takes into account momentum conservation but the average kinetic energy grows to 
infinity for t — > 00 and hence the process of ultimate relaxation to equilibrium is not properly described. However, it is 
expected that in the limit of large mass M and for slowly varying potential V^(x) the decoherence time is much shorter 
than the energy dissipation time scale. Therefore the master equation (35) can be considered as a good approximation 
for the study of pure decoherence in the relevant regime of large body at slowly varying potential and rare collisions 
(see Section III. D.3). 
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4- Bloch-Boltzmann equation 



In quantum optics the investigation of an active atom interacting with radiation and immersed in the thermal bath 
of perturbers (typically , noble gas atoms) is a very important topic [18]. To simplify the description one assumes 
that the spatially uniform distribution of active atoms and strong decoherence due to collisions with perturbers justify 
the use of the density operators diagonal in momentum (or velocity) representation. Therefore, we consider a family 
p(v) > of velocity-dependent positively defined operators (matrices) acting on the Hilbert space describing the 
internal structure of an active n-level atom. Using a standard form of the master equation (13) one can easily derive 
the most general evolution equation for QDS which preserves such partially diagonal structure of the density matrix 
[19]. It reads 




where {S a } form a linear basis in the space ofnxn matrices and 

foranyv,v', /C a6 (v;v') is positively defined , 7 a &(v) = / rf 3 v'/C ha (v; v') . (37) 

Equation (36) provides an example of a classical-quantum description which combines a generalisation to n-level case 
of the MME (22), called Bloch equation with a classical linear Boltzmann equation for spatially uniform case. Similar 
type of equations found applications in the quantum measurement theory [20] . 

D. Ito-Schrodinger equations 

In the classical theory the Fokker-Planck equation describing Brownian motion in terms of probability distribution 
over the phase-space of Brownian particle is completely equivalent to the Langevin equation (in Ornstein-Uhlenbeck 
formulation) for the particle's random trajectory. The later equation is a prototype of stochastic differential equation 
which can be written either in Ito or Stratonovitch form [1]. The similar equivalent description exists for quantum 
dynamical semigroups with the generator (13) in terms of the following Ito-Schrodinger equation 

#(i) = -iHip{t)dt V o* V J^)dt - i Vj%l>(t)dBj (t) . (38) 

3 3 

In (38) Bj(t) denote independent and normalized Brownian motion processes and heuristically "4f" corresponds 
to "white noise". The short-hand notation df(t) = a(t)dt + b{t)dB{t) has a meaning of the corresponding integral 
equation 

f(ti) - f(t 2 ) = f 2 a(t)dt + f 2 b(t)dB(t) (39) 

where the second integral is an " nonanticipating" Ito integral given by the probabilistic limit of the sums 
^2 b(sk)[B(sk+i) — B(sk)] with t\ = si < s 2 < ■ ■ ■ < s n = i 2 - To perform computations using stochastic inte- 
grals one applies Ito rule 

dB i {t)dB j {t)=6 ij dt (40) 

and E[dBj(t)] = where E[-\ denotes the average with respect to Brownian motions. One should remember that 
dBj (t) — Bj (t + dt) — Bj (t) is always statistically independent on the solution of stochastic differential equation up to 
time t. Applying those rules to the density matrix of the open system defined as an average p(t) = E[\ip(t) >< ip(t)\] 
we obtain 

d Pt - dE[\ip(t) >< V(i)l] = E[\d*l>(t) >< Ht)\ + IV>W >< #(*)| + |#(i) >< #(t)|] 

= -i[H, Pt ]dt - i{]T V*V„ Pt }dt + V jPt V*dt (41) 

3 3 
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that coincides with MME (13). The Ito-Schrodinger equation (38) for a general case was firstly proposed in [21,22]and 
then derived by different authors as a useful tool for numerical simulations of quantum open systems [23] . For such 
applications it is convenient to use an equivalent nonlinear modification of (38) with the normalized solutions <f)(t) = 
i/j(t)/\\ip(t)\\. One should notice that to obtain (41) we can replace the Brownian motions Bj(t) by arbitrary stochastic 
processes Xj(t) ( continuous or jump processes) with independent increments such that E[dXi(t)dX j(t)] = Sijdt. The 
different processes represent unitarily equivalent representations of the singular reservoirs constructed in terms of 
"time ordered Fock spaces" [21]. Possible physical interpretations of these different representations are discussed in 
[24]. 

An alternative approach in terms of "quantum noises" has been proposed in [25] and subsequently developed into 
a mathematically rigorous theory of noncommutative stochastic differential equations. 

E. Problems and pitfalls 

Although, the axiomatic approach allows to write down the MME (13,14) in a mathematically consistent form, we 
generally need additional phenomcnological insights which could lead to concrete form of the operators {Vj}. The 
first problem is the relation between the dissipative part in (13,14) and the Hamiltonian of the open system. Here 
the demanded structure of the stationary state and detailed balance condition arc helpful [7]. For example adding a 
nonlinear term to the Hamiltonian of the harmonic oscillator we have to modify the dissipative part too, in contrast 
to often used simplified models with linear dissipation. The symmetries of the system and reservoir give another 
restrictions on the form of generators [26] . 

For many interesting physical problems we have to include time dependent-external fields acting on the system. 
Very often it is done by simply adding the proper time-dependent term to the Hamiltonian without changing the 
dissipative part. This can be justified, however, for weak and quickly oscillating fields while in general case the whole 
Markovian approximation breaks down [27] . Another extreme case is a slowly varying external potential which leads 
to a proper time-dependent generator which for any time reproduces, for example, return to an equilibrium state 
given by a total instant Hamiltonian [28] . 

A particular care is needed when we put two open systems Si, S2 with their "privat" reservoirs Ri, R2 into contact 
by adding their interaction Hamiltonian V12. In this case we cannot simply add two dissipative generators obtained 
previously for decoupled systens 5*1,52. This is obvious for the case of two heath baths at the same temperature 
which should drive the interacting system 5i + 52 into its Gibbs state which of course depends on the new element 
Vi 2 . 

III. CONSTRUCTIVE APPROACH 

The goal of the constructive approach is to derive, using mathematically controlled approximations, Markovian 
dynamics of the quantum system interacting with the quantum environment. We restrict ourselves to the deriva- 
tions which respect complete positivity of the reduced dynamics. In author's opinion the very scheme of quantum 
open systems approach and particularly Markovian approximations are consistent only under well defined conditions 
discussed below. 

A. Basic assumptions 

In the theory of quantum open systems we decompose the Universe into three parts (5, 0,R). The first one is an 
open system 5 controlled by an "active" observer O and interacting with the "passive" rest R which is called reservoir. 
If the state of the system 5 + R is given by the density matrix psr that we asume that the reduced density matrix 
ps = T^trPsr possesses an operational meaning i.e. that the observer O can prepare an arbitrary initial state of 5 at 
a certain moment i = and than measure the evolved state at the moment t > 0. We often assume also that the 
observer can control to some extend the evolution of 5 switching-on external time-dependent " fields" . The physical 
structure of O and the description of measurements performed na 5 is a subject of quantum measurement theory and 
will be not discussed here. 

One should notice that the above scheme contains a tacit assumption that the interaction between 5 and R is 
weak. Any preparation of the initial state of 5 or its state measurement takes time t p = H/5E where SE is an 
energy resolution of state preparation or measurement. For a quantum system 5 with a discrete spectrum SE should 
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be smaller that the typical separation between the energy levels while for systems with continuous spectrum (e.g. 
quantum Brownian particle) t p is directly related to a chosen time-scale of observation (see Section III.D.3). Obviously, 
during preparation or measurement processes the system S should not be strongly perturbed by the environment R 
what implies the inequality 

t P «T S (42) 

where ts is a characteristic time scale of relaxation of S due to the interaction with R. 

The first step towards mathematical description of the open system consists in writing down the formal "bare" 
Hamiltonian 

H SR = H s + \Y,S k ®R° k +H R (43) 

k 

where Hg is a bare Hamiltonian of S, the interaction term is proportional to a dimensionless coupling constant A and 
contains bare reservoir's operators For some applications we shall add a time-dependent contribution V(t) to the 
Hamiltonian Hg in order to take into account external fields used by the observer O to control our system S. 
The next important assumption is the product structure of the initial state for S + R 

Psr(0) = P®u r (44) 

with an arbitrarily prepared initial state of 5* and a fixed reference state of R. Again this assumption can be justified 
if the interaction is weak in a defined below sence [29] . Namely, the system S + R for times preceding the preparation 
at i — is expected to be at equlibrium or ground state wsr- In order to produce the initial state (44) by applying 
external perturbations to S we must assume that 

u S r ^ ui s ® u) R . (45) 

The state of S for times t >> t p is then given by the reduced density matrix (we omit subscript S) 

pt = Tr R (U tP ® w R U*) = A t (p) (46) 

where U t = cxp{—itHg R } and Tr^ is a partial trace over the Hilbert space of R. The dynamical maps {A t ;t > 0} 
are by construction completely positive and trace preserving. 

However, the practical implementation of the above prescription leads to difficulties. Take as an example, a model of 
atom interacting with electromagnetic field at zero temperature. The corresponding Hamiltonian (43) is well-defined 
only if an ultraviolet cut-off at the frequency uj c is introduced replacing "bare" operators by cut-off dependent 
regularized Rk- Then a simple lowest order evaluation of (46) shows that the interaction produces a cut-off dependend 
Hamiltonian correction to the bare Hamiltonian Hg. When lo c — > oo corrections diverge and must be cancelled by 
proper counterterms in the bare Hamiltonian (renormalization procedure). In the case of our example it leads to mass 
and charge renormalization and slight modification of the Coulomb potential which influences energy levels of the 
atom (Lamb-shift) [30]. A different example is an electron immersed into a polar solid. The interaction with lattice 
ions modelled by the elect ron-phonon coupling (with natural well-defined cut-off) completely changes the propeties 
of an electron leading to the so-called polaron model. In both cases the initial bare system is strongly modified by a 
"cloud" which consists of certain degrees of freedom of the reservoir [31]. Hence the decomposition into S + R is not 
unique and the practical criterion should be the weakness of the residual interaction between the " dressed" system 5" 
and the rest of the reservoir R'. Practically, we describe the sytems S', R' by the same Hilbert spaces as S, R and the 
Hamiltonian of the bath H R remains the same too. We introduce, however, a physical (renormalizcd) Hamiltonian Hs 
and modify the structure of interaction Hamiltonian introducing suitable formfactors with cut-offs and often leaving 
only resonant terms ("rotating wave approximation" in quantum optics, Lee models in particle physics) to reduce 
creation of the "cloud". We shall ilustrate these problems later on in the case of 2-level open system (see Sections 
III.C.2, III.D.2). 

B. Models of reservoirs 

Similarly to the classical case represented by the standard example of Brownian motion the details of the reservoir 
should not be essential for the dynamical behaviour of an open system S at least on the proper "coarse-grained" 
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time scale. Instead we expect that the mechanisms illustrated by limit theorems in classical probability should work 
in the quantum domain also. Indeed, take a sequence of quantum observables A\, A2, An which are statistically 
independent (i.e. < A k Ai >=< Ak >< Ai >) or even weakly dependent. Then, one can prove a quantum version of 
the Central Limit Theorem for TV — > 00 

1 N 

-j= J2(A 3 - < Aj >) a(f A ) + a+(f A ) (47) 

where a(./U), o- + {Ja) are "smeared" bosonic anihilation and creation operators [32,33]. The limit is understand in the 
sense of convergence of all correlation functions, where for Bose fields we choose a "vacuum" as a reference state. 
The detailed structure of the Fock space for bosonic field is discussed in [10,32,33] and is not relevant here. One 
should think about bosonic fields as quantum counterparts of classical Gaussian random fields and the " vacuum" can 
represent any quasi- free state of bosons as for example arbitrary Gibbs state for noninteracting Bose gas. The index 
j for Aj can correspond both to different parts of the reservoir and to different time instants . The limit theorem (47) 
reflects the physical picture of influence from R on S which can be decomposed into a sum of many almost independent 
small contributions. Such case is usually called diffusive regime and due to quantum central limit theorem (47) is well 
described by the harmonic oscillators bath with a weak linear coupling in the interaction Hamiltonian 

HTr = A £ Sfc ® HZ*) + ft+ (A)] • ( 48 ) 
k 

Diffusive regime does not describe all physically interesting situations. Another, low density regime corresponds to 
rare events of scattering by essentially independent "particles" from R each of the scattering processes takes short time 
but need not to be weak. This is similar to classical limit theorems leading to Poisson distributions and processses. 
The physical model of such reservoir is a free fcrmionic or bosonic gas and the bilinear interaction given by 

HTr = E S * ® ia + (.fk)a(g k ) + a+(g k )a(f k )} (49) 

k 

with a small parameter being now the density of gas ^< a + (fk)a(fk) >■ 

For some applications the quantum nature of the reservoir is not very important, like for instance in the high 
temperature regime. In this case one can replace a quantum reservoir R by a time dependent random Hamiltonian or 
in other words operator-valued (self-adjoit) stochastic process. A particular choice of white-noise process leads to a 
stochastic Ito-Schrodinger equation (38) with Vj = . 

C. Markovian limits 

The general reduced dynamics (46) does not satisfy the semigroup composition law (7). However, one expects 
that at least for a certain coarse-grained time scale, roughly determined by t$ (see (42)), this law is often a very 
good approximation. Physically, it is true if the exact state of the system S + R given by U t p ® ujrU^ does not 
differ locally from the state A t (p) ® u)r. Here "locality" is determined by the radius of interaction between S and R. 
Mathematically this condition can be expressed as a sufficiently fast decay of the reservoirs correlation functions 

R u {t) = Tv{u R R k (s + t)ifc(a)) (50) 

where Rk{t) = exp(itH r) Rk cxp(— it Hr) and we assume that the reservoir's reference state ojr is stationary with 
respect to its evolution. We shall argue that those correlation functions contain the total information about reservoirs 
which is relevant for different Markovian regimes. Introducing the reservoir's relaxation time scale tr we may write 
the standard condition for the validity of the Markovian approximation as 

tr « t s . (51) 

However, as we shall see, important examples of reservoirs may not have a natural decay time scale tr and one has 
to include the averaging effect of the Hamiltonian dynamics of S in order to satisfy (51). Moreover, rigorous analysis 
shows that the Markovian approximation needs also sufficcntly rapid decay of higher order multitime correlation 
functions which can be proved for quasi-free models of reservoirs discussed in the previous subsection. 
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In the present review we are not going to reproduce involved rigorous derivations which can be found in the literature 
[34-36]. Our aim is to discuss physical assumptions behind different Markovian regimes and proper approximation 
schemes based on the relevant order of perturbation which yield mathematically consistent results. For notational 
simplicity we use as an illustration the model of an effective 2-level open system with a physical (renormalizcd) 
Hamiltonian H$ and the interaction one Hi nt of the form 

H s = ^ecr 3 , H int = Acti <g> R (52) 
where Ofe, k = 1,2,3 are Pauli matrices. We always assume that 

Tt(lo r R) = . (53) 
In the heuristic derivations of MME the lowest order expansion for quantum dynamical semigroup of the form 

Pt = e- itH * (P + E /' - \{V*{s)V j {s),p}ds)e itH ^ + 0(t 2 ) (54) 

is compared with the relevant term in the expansion for the reduced dynamics (46). More precisely, it is enough to 
compare only the transition map term (" gain" ) 

P^Hi Vj(s)pV*(s)ds (55) 

3 J ° 

with the corresponding expression obtained from (46) using different assumption concerning mainly important time 
scales in the joint system S + R. Any Hamiltonian-type corrections are put equal to zero according to the renormal- 
ization procedure. 

1. Singular coupling or white noise Anzatz 

Using second order Dyson expansion for the total dynamics of S + R with the interaction Hamiltonian treated as 
a perturbation we obtain the following formula for the manifestly completely positive "transition map" which has to 
be compared with (55) 

p^> I dsi / ds 2 r Tr R (H ln t{si)p®uj R H m t{s2)) ■ (56) 
Jo Jo 

For our 2-level system we obtain 
P^ X2 J o ds iJ o rfs 2 i?(s 2 - Sl )(e ie ( S2 - Sl V-pa + + e- ie ( S2 - Sl V + p CT -+e- te(s2+Sl) cr-p C r-+e i£ ( S2+Sl V+pa+) . (57) 

The simplest approximation which produces the Markovian transition map term (55) is based on the assumption that 
the correlation time of the reservoir is much shorter than any other relevant time scale i.e. 

t r « min{T S ,T H } (58) 

where tr denotes Heisenberg time scale for S related to the average energy level spacing AE by th = (AE)^ 1 . Under 
the condition (58) we can use the white-noise Anzatz 

/oo 
R{t)e Uul dt (59) 
-OO 

which leads to the QDS-type of the transition map corresponding to the following MME 

j t P = -\^,p]-—R{Q)^u[<yup\] ■ (60) 
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For a general case the substitution of the type (59) leads to MME (13,14) with Vj = V* which can be rewritten in a 
double commutator form. Moreover, the dissipative part does not depend on the Hamiltonian one and we can even 
add a time-dependent term V(t) to obtain 

jtPt^^KHs + vmpt}-^^'^^^ • ( 61 ) 

3 

The dynamics generated by (61) is very special because the dynamical maps A t are bistochastic, i.e. A t l = 1. It 
means, depending on the physical context, that the infinite temperature or microcanonical state is preserved and the 
H-theorem holds i.e. for any solution of (61) [37] 

S(p t ) > S(p s ) ,for t > s > , S(p) = -Trplnp . (62) 

There are several instances where the equations of type (61) are useful and provide reasonable approximations to the 
exact dynamics of the open system. These are the cases where the Hamiltonian self-evolution due to H$ + V(t) is very 
slow or H$ + V(t) commutes with Sk in (43)and therefore the pure decoherence processes dominate over the energy 
redistribution. Another case is the high temperature limit where T » E max and E max is a maximal difference of 
energy levels of the system. The white noise Anzatz can be obtain from the Hamiltonian models by limit procedure 
called singular coupling limit which leads to unphysical limiting reservoirs with unbounded from below Hamiltonians 
[35]. 

One should remember, however, that the lack of memory expressed by the (5-like correlations in (59) contradicts the 
quantum nature of the reservoir (see Section III.D.l) and the noise governing the dissipative part of (61) is essentially 
classical. Hence, in author's opinion, equations of the type (61) are too rough to describe properly the control of 
decoherence and dissipation in the context of quantum computations and error correction schemes. The same is true 
for their discrete-time versions which often appear in the literature on this topic [9] . 



2. Weak coupling limit 



Another approximation scheme takes into account the interplay between the Hamiltonian dynamics of S governed 
by H$ and the coupling to R. It allows to describe, for instance, the equilibration process leading in the case of a 
heat bath to a final Gibbs state of S. We assume the following relations between the relevant time scales 

TH « T S , , T R « T S ■ (63) 

As a consequence of (63) in the example (57) the last two non-resonant terms can be omitted and the obtained 
transition map corresponds to the following MME 

j fPt = - l -e[a 3 ,p t ] ^^R(e)([a-,p t a+] + [a~ p t ,a+}) + R(-e)([a+ , p t a~] + [aW"])) . (64) 

which is a special case of (22) with 5 = 0. The similar procedure applied to a harmonic oscilator linearly coupled to 
an environment yields MME (27) 

For a general case (43) the MME obtained under the conditions (63) reads 

j t p t = -i[H s ,p t ] + y E Rki{u)([Sk(w),(HSt(u>)] + [S k (w)(H,St(w)]) (65) 

uj-.k.l 

where 

S k (t) = e itHs S k e- itHs = £ SkMe-- 4 (66) 

UJ 

and 

/oo 
R kl {t)e- itu dt . (67) 
-oo 

Obviously, as for any eigenfrequency lo the matrix \Rki(u))] is positively defined we can always rewrite (65) in a 
diagonal form (13,14). 
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The approximation procedure leading to MME (65) ( in the interaction picture) can be made mathematically 
rigorous introducing the concept of weak coupling or van Hove limit applied to the reduced dynamics in the interaction 
picture which consists of using the rescaled time r = \ 2 t for A — > 0. The details can be found in [34] see also the 
related idea of stochastic limit in [33] . 

MME of the type (65) possess several interesting properties. The dissipative part of the generator commutes with 
the Hamiltonian one and the diagonal elements of p t in the energy representation evolve independently of the off- 
diagonal ones. The corresponding transition probabilities which appear in the Pauli master equation for diagonal 
elements are exactly equal to those calculated using the Fermi Golden Rule [38]. In the frequently used case of the 
reservoir being at thermal equilibrium (heat bath) the KMS relation of the form 

Rki(-iu) = e-^Riuiw) , (3 = T- 1 (68) 

is valid and implies that that the Gibbs state p eq = Z^ 1 cxp{— (iH$} is an invariant state for MME (65). If for all 
k, u, [Sk{oj), X] = implies X = cl the Gibbs state is ergodic i.e. for t — ► oo and any initial state the solution of (65) 
Pt — * Peq- Moreover, the Pauli master equation for the diagonal elements of p t satisfies detailed balance condition. 
The part of the generator in (65) which is given by terms with u> = does not influence the diagonal elements of p t 
and describes pure decoherence which is not accompanied by the energy change. 

One should notice the important differences between the regimes characterized by inequalities (58) and (63) re- 
spectively. The former allows white-noise Anzatz and hence the strictly Markovian approximation entirely due to the 
memoryless reservoir. In the later case Markovian approximation is valid for the time evolution averaged over many 
periods of Hamiltonian evolution of the open system. It follows that the dissipative part of the semigroup generator 
strongly depends of the Hamiltonian Hs and therefore we cannot freely add a time dependent part V(t) to H$. 



3. Weak coupling and thermodynamics 

We have already seen that the MME (65) obtained using weak coupling limit satisfies the 0-law of thermodynamics, 
namely the heat bath drives an open system S to its equlibrium state at the same temperature. Introducing a 
time dependent Hamiltonian Hs(t) which varies on the time-scale much longer than all other time-scales in (63) we 
can apply the same method used to derive MME (65). As a result of such combined weak-coupling and adiabatic 
approximation we obtain an inhomogencous in time MME of the form 

j t Pt = -i[H s (t),p t ] + L D (t) Pt = L{t) Pt , t > . (69) 

where the dissipative part L D {t) is computed according to the formulas (65,67) with the Hamiltonian Hs replaced by 
its temporal value Hs(t). We can now formulate the 1st law of thermodynamics as 

where E(t) = Ti(p t Hs(t)) is an internal energy of S, W(t) = f*Tr(p s ±H s (s))ds is the work performed on S by 

external forces and Q(t) = J Q Tt [(-^p s )Hs(s)]ds is the heat supplied to S by R [28,39]. 

In order to illustrate the second law of thermodynamics we use the relative entropy for a pair of density matrices 
p,a 

S(p\<j) =Tr(plnp-plno-) (71) 
and the following inequality valid for any trace preserving completely positive map A [37] 

S(p\o-) > S(Ap\Aa) . (72) 
For the time evolution governed by the generator L(t) (eq.(69)) which posseses a temporal stationary state 

L(t) Peq (t) = , Peq (t) = Z- 1 c^{-(3H s (t)} (73) 
the inequality (72) implies the following form of the second law of thermodynamics for open systems 
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where o~[pt] > is an entropy production and the second term describes the entropy exchande with the heat baths. 
The formula (74) can be easily generalized to the case of open system coupled to several heath baths at different 
temperatures [28,40]. 

The presented derivation of the three fundamental laws of thermodynamics has been possible within the assumptions 
of weak coupling between S to R and the adiabatic change of external forces. Beyond this regime, even the unique 
definitions of the fundamental thermodynamic notions are not obvious. One should mention that the deviations from 
the thermodynamical behaviour for strongly coupled quantum systems attracted, recently, attention of the scientific 
community [41]. 



4- Low density limit 

Another situation which justifies the use of MME is an open system S with the discrete spectrum Hamiltonian 
Hs = X)fc e k\k >< k\ interacting with a dilute gas. The rigorous results can be found in [33,36] while here we present 
a heuristic approach. Consider first an abstract formulation of the impact approximation for the scattering problem 
given by the Hamiltonian H = H + V. Starting with the identity 

U t = e- ltH = e - UHo (l + e tsHo Ve- tsH ds^ = e- ltH °(l + e lsH °V[e- lsH e lsHo ] e - lsH ° ds"j (75) 

we replace e - lsH e lsH o by fi + = lim s ^ 00 e ~ tsH e lsH " what makes sense for t » t co u where t co u is a typical collision 
time for the discussed model. Hence we have 

U t « e~ ltHo (} + J e lsHo Vn + e- lsHa ds^ (76) 

and instead of the transition map (56) we obtain its low density counterpart 

p~ [ d*i / ds 2 Ti R (T( Sl )p(g,uj R T(s 2 )) . (77) 
Jo Jo 

where T(s) = e tsH oTe~ lsH ° and T = VSl + . In the next step we obtain a transition map similar to that derived in the 
weak coupling limit by averaging over the Hamiltonian evolution of S and eliminating the non-resonant oscillating 
terms. Take a single particle of the bath being in a state described by the density matrix 



ujr = r 3 f d 3 P G(p)|p X P | , < p|p' >= S 3 (p - p') 

JR3 



(78) 



normalized in a cube £ 3 and diagonal in the momentum representation. Putting H = Hs + J R3 d 3 p E p \p >< p| 
we obtain from (77) and for t >> t co u , th the final form of the transition map which commutes with Hamiltonian 
evolution 

p^tr 3 Y^( d 3 p [ d 3 pG{p)^5{E p ,-E p +uj)T UJ {p,p')pT:{p,p') (79) 
^7 iR3 Jr 3 

where 

2L(P,P')= E <k,p'\T\l,p>\k><l\ . (80) 



The transition map (79) corresponds to a single-particle scattering during the time interval t. Because in the volume 
£ 3 we have N particles the low density approximation implies additive effect producing the factor v = N/t 3 (instead 
of £~ 3 ) which remains finite in the thermodynamical limit. Therefore, the final form of the MME valid under the 
conditions 



th « t s , ,t coU « t s . (81) 

is the following 
d 

p t = ~i[H s ,pt\+vn- } j 

/r 6 



dt 



Pt = -i[H s , P t] +™Y,J d3pd3p ' G ^ S ( E p' - e p + ") (V-(P> P'). PtT^P, P')] + PL(P, P'K *Z(P, P')]) (82) 
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Similarly to the weak coupling limit the rigorous derivation of MME (82) involves the limit procedure v — > for the 
interaction picture version of the reduced dynamics with the rescaled time r = vt [36]. 

It is not difficult to show that the MME (82) has similar properties to MME (65). Namely, the diagonal elements 
of pt evolve independently of the off-diagonal ones, for the equilibrium momentum distribution of scatterers G(p) ~ 
exp(— /3E p ) the Gibbs state p eq ~ exp(— (3Hs) is stationary and under natural conditions ergodic. The discussion of 
the Section III.C.3 can be applied also for inhomogeneous in time versions of MME (82). 

5. Application to Block- Boltzmann equation 

The MME (82) can be a starting point for the heuristic derivation of the Bloch-Boltzmann equation (36) [19]. The 
active atom with finite number of internal levels can be put into a finite box to take an advantage of the dicrete 
spectrum of its momentum and kinetic energy operators. To such a discretized system of an active atom we can apply 
the derivation of the previous section. As an initial state we choose a quasi-diagonal density matrix 

p = E^ v )i v><v i ( 83 ) 

V 

with respect to the orthonormal velocity basis {|v >}. Under some reasonable conditions the quasi-diagonal density 
matrices remains quasi-diagonal during the time evolution and finally we can go with the size of the box to infinite to 
obtain continuous velocity spectrum. The final result of this procedure which involves also the transition to the center 
of motion reference frame is the Bloch-Boltzmann equation (36) (37) with (we put here fi for readers convenience) 

n 

H s = Y, e M >< j\ ' \- H s> S ^ = ' a = 2 ' -' n2 ( 84 ) 
/C a6 (v; v') = (2n) i h 2 p- 3 N p S LJa . LJb / dV r [ d 3 v r 5 3 {v - v' - ^(v r - v' r ))s(^(v 2 r - v' r 2 ) + hcj a ) 



xTU(v'-v;)T Q (v r ,v;)T h (v r ,v;) (85) 

where v r ,v£. denote relative velocities, W p (v) is the equilibrium velocity distributions of perturbers, to is a mass 
of an active atom, /i is a reduced mass of an "active atom - pcrturbcr" system, N p is a density of perturbers and 
the functions T a (., .) are related to the T-matrix calculated in the center of motion reference frame by the following 
expression 

n 

^T a (v r ,v;)5 a = E < VrJ\T\v' r ,f > \j >< j'\ ■ (86) 

a 3,j'=l 

Due to the presence of Kronecker's delta 5 UatUb in (85) the matrix /C a b(v;v') is positively defined as demanded by 
the conditions (37). One should mention that the obtained Bloch-Boltzmann differs from the existing ones which 
generally do not preserve positivity of p(v) [18]. 

D. Problems and pitfalls 

The derivations of the MME presented in the previous sections are based on certain assumption concerning the 
decay of correlations in the environment, separation of energy levels of the open systems and the magnitude of the 
coupling constant or the density of perturbers. These conditions are formulated in terms of relations between the 
different time scales (42), (51), (63), (81). Now, we discuss briefly some of them and refer to the important physical 
situations beyond the presented scheme stressing the related difficulties and misconceptions. 
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1. Memory effects 



The most demanding condition for the lack of memory of the reservoir is an inequality (58) which allows the white 
noise Anzatz Rki(t) ~ Rki(0)5(t). It means that the spectral density matrix Rki(u>) is weakly dependent on u> in the 
relevant energy region. Such an assumption essentially contradicts the KMS condition (68) what is the source of the 
so-called thermal quantum memory characterized by the time scale tt = T _1 . 

Another difficulty with the assumption of memoryless reservoir can be illustrated by a model of 2-level atom coupled 
to the quantum electromagnetic field at the vacuum state (T = 0). Puting in cq.(52) the standard "dipolx electric 
field" interaction one obtains 

~ 7 nZ > ~ u^e'"'"' (87) 



(t + iuo, 



where ui c is a cut-off frequency assumed to be larger than any energy scale relevant for this model. The decay of 
reservoir's correlations is powerlike and does not possess any natural time scale. Moreover, the correlation function 
is singular at the origin and i?(0) = for the removed cut-off (u> c — > oo). Fortunately, in the weak coupling regime, 
the Markovian behaviour can be restored on the coarse-grained time scale determined by ts - spontaneous emission 
time. This is due to the averaging effect of fast Hamiltonian oscillations which allow to replace R(t) by R(t)e let - the 
function which effectively acts like R(e)5(t) under the integral in (57). The situation is different when Hs is replaced 
by a time-dependent Hamiltonian and /or collective effect for multi-atomic systems are relevant. Namely, that fast 
variations of the Hamiltonian introduce high frequency contributions increasing the decay rates (see eq.(87)) and even 
a system of two atoms possesses degenerated energy levels for which averaging effect described above does not apply. 
All that implies serious limitations on the use of Markovian approximation in the context of controled quantum open 
systems, the problem which is crucial for quantum information processing [42]. 



2. Decoherence vs. dissipation 

The dynamics of an isolated quantum system is governed by its Hamiltonian Hs and is characterized by two 
fundamental features : initial pure states remain pure , the average energy is a constant of motion. On the contrary, 
for an open quantum system S interacting with a quantum environment R which starts its joint evolution from the 
product state, the entangled states of S+R are developed in the course of time what lead to the appereance of reduced 
mixed states of S and the energy exchange between S and R. The first phenomenon is called decoherence and the 
second one dissipation. Although decoherence and dissipation arc usually present at the same time, model calculations 
show that for large quantum systems approaching the border between quantum and classical worlds decoherence acts 
on a much faster time scale than dissipation [43,44,17]. As decoherence seems to be a more important agent in the 
context of quantum measurement theory and quantum information processing [9] it is convenient to discuss models 
describing pure decoherence often called dephasing which is not accompanied by the energy exchange. 

The pure decoherence is described by the models with Hamiltonians (43) satisfying the condition [H s , Sk] — what 
implies the same condition for the renormalized Hamiltonian [Hs,Sk] — for all k. In Markovian approximations 
the decoherence rates are proportional to Rki{0) (see (59,60) and (65-67)) which is typically zero for the systems 
linearly coupled to bosonic fields . This is the case for electromagnetic interaction (87), the same result holds for 
linear coupling to phonons. In a general case one can formulate the following "no-go theorem" for pure decoherence: 

For quantum open systems linearly coupled to bosonic reservoirs decoherence is always accompanied by dissipation. 

The "physical proof is rather simple. Any irreversible decoherence must be related to an irreversible change in 
the environment. In the case of linear coupling to bosonic field (48) this change can be realized only by emission 
or absorption of a boson - the process which changes energy of an open system as well. Obviously, in the case of 
scattering process governed by the bilinear interaction Hamiltonian (49) we can alter other quantum numbers of the 
environment's state (e.g. its momentum) keeping the energy of S conserved (elastic scattering). 

One can easily find in the literature the models of dephasing based on the linear coupling to bosonic field which are 
essentially variations of the so-called Caldeira-Leggett model [45-47] . To explain this apparent contradiction with our 
"no-go theorem" consider the simplest version of the spin - boson model defined by the Hamiltonian (here Hs = 0) 

POO poo 

H S R = Xo- 3 (g> dw[f(oj)a(u) + f(u)a + (uj)}+ / doj lu a + (uj)a(uj) (88) 
Jo Jo 
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acting on the Hilbert space 



Hsr = C 2 ® T B (L 2 [0, oo)) ) = T B (L 2 [0, oo)) ) © .Fb (L 2 [0, oo)) (89) 

where .Fb(.L 2 [0, oo))) is a bosonic Fock space over a single-particle Hilbert space L 2 [0,oo) and [a(u)),a + (w')] = 
6(uj — u'). The unitary Weyl operator U g acting on fields operators as 

U g a(u)U; = a(u)+g(w) (90) 

exists if and only if g e L 2 [0, oo). Puting 

g(u) = \Lo-\f{iu) (91) 

we can diagonalize Hsr (88) 



UgH S R\J* g = I duo uj a + (uj)a{uj) + const. (92) 

where 



o 



The ground states subspace of the diagonalizcd Hamiltonian (92) is spanned by the vectors |1 > ®|0 >, |2 > > 
and therefore the corresponding degenerated ground states of Hsr are given by 

|1 > ®\4>[-g] > , |1 > <8#[ff] > , where >= C^fi > , g(u>) = Aw" 1 /^) . (94) 

The vectors |</>[±g] > are coherent states in Tb (L 2 [0, oo)) and 

<0[- 3 ],4 9 ]>=cxp{-2||. 9 || 2 } . (95) 

Assume now that we would like to describe dephasing of our spin using the model given by (88). Then we should have 
R(0) — A 2 |/(0)| 2 > what implies due to (91) that the function g(u) is not square integrable i.e. \\g\\ = oo. The same 
divergence appears for any "ohmic" or "subohmic" coupling |/(cj)| 2 ~ cu s around with < s < 1. It follows that 
the diagonalizing transformation of a bosonic field (90) (91) cannot be implemented by the unitary operator on the 
Fock space and therefore the formal expression (88) does not define a meaningful bounded from above Hamiltonian. 
It means that for an arbitrary coupling constant A the model given by (88) is either nonphysical or cannot describe 
exponential (Markovian) dephasing. In a less mathematical language the large value of \\g\\ due either to a large 
coupling constant or to a large integral J °° dui o>~ 2 |/(cj)| 2 makes the coherent states |0[±g] > corresponding to a 
"cloud" and the vacuum |f2 > almost orthogonal. As a consequence the standard choice of the initial state as a 
product state p® \tt >< Q\ is inappropriate and imposiblc to prepare. The proper initial state should have a support 
spanned by the "dressed" ground states (94). The computed lost of coherence for the former choice of the initial state 
is therefore unphysical and describes the spurious process of a " cloud formation" . Again the problem of a proper 
decomposition of the total system into open system and environment such that the effective interaction between them 
is weak, is crucial for the physical interpretation of the obtained results. 

One should mention that from the mathematical point of view the above example illustrates the subtle problem 
of nonequivalent representations of canonical commutation relations for systems with infinite number of degrees of 
freedom [48]. 



3. Open systems with continuous spectrum 

In the derivations of MME based on the weak coupling or low density limits the discretness of Hs spectrum plays a 
crucial role. The averaging over Hamiltonian oscilations justifies the Markovian approximation on the coarse-grained 
time scale and the canceling of non-resonant terms is a necessary condition to preserve complete positivity of the 
QDS. This interplay between the self-evolution of S and the interaction with R produces the desired properties of 
QDS like for instance the relaxation to a proper equilibrium state for R being a heat bath. 
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However, there are important examples of open systems with continuous spectrum of Hs, the most studied is a 
quantum Brownian particle. For simplicity, we discuss first the case of Brownian motion in one-dimensional space. 
The most frequently used MME for this case is the so-called Caldeira-Leggett equation [47] of the form 

j t Pt = -i[H s , Pt] - h[X, {P, Pt }} - 2M 1 T[X, [X, p t ]] (96) 

where X, P are position and momentum operators, Hs — P 2 /2M + V(X), 7 is a friction constant and M- mass of the 
Brownian particle. V(X) is a generic potential which can produce both continuous and discrete parts of the energy 
spectrum. 

The MME (96) possesses the following well-known drawbacks: 

1) the solution of (96) does not preserve positivity of the density matrix, 

2) the Gibbs state (perhaps unnormalized) ~ exp{— Hs/T} is not its stationary state. 

The first drawback can be cured by adding the term — 7(8MT)~ 1 [P, [P, p t ]] which allows to write the corrected 
MME in a standard form (13,14). However the dissipative part of this new generator corresponds to a damped 
harmonic oscillator one (27) (with 71 = 0) and the Hamiltonian part gains the correction ~ (X P + PX) with a rather 
unclear physical interpretation. 

In author's opinion the above difficulties have its source in the underlying Hamiltonian model for S + R system, 
the so-called Caldeira-Leggett model, which is essentially a non-zero temperature version of the model (88) with 03 
replaced by the operator X and the "forbidden" ohmic choice of |/(w)| 2 ~ u> around zero. In other words according to 
the "no-go theorem" the equation (96) which in the extreme heavy particle limit (M — > oo)describes pure decoherence 
cannot be derived from a physically admissible model with linear coupling to the bosonic reservoir. 

It seems, that for a continuous spectrum of Hs the lack of a natural time scale provided by the Hamiltonian 
evolution makes impossible to find a single MME which accurately describes all relevant stages of the evolution of S 
[49]. We do not mean here the well-known and expected deviations from the Markovian (exponential) behaviour for 
very short times typically ~ lo^ 1 (ui c - cut-off frequency) and very long ones (due to the boundness from below of 
the Hamiltonian of R). In the simplest case of a Brownian particle in the free space described by the center of mass 
position X and the total momentum P the situation can be summarized in a following way. 

The proper model of an environment should involve interactions of the type (49) describing scattering process with 
particles of the medium. Beside the direct collisions with atoms, molecules, photons etc. the other "bilinear" processes 
are possible within this model like for instance absorption of a foton followed by the excitation of the internal degrees 
of freedom of the Brownian particle and the time-reversed process [17]. 

The following different but approximatively Markovian stages of the evolution can be singled out: 

1) Pure decoherence stage, when the decay of macroscopically distinguishable quantum superpositions into mixed 
states dominates over energy thermalization. For the low density medium and/or small Brownian particle this stage 
is well-described by the eq. (35) on the time scale determined by Tdec — (/ d 3 kn(k)) . For dense media or/and 
large particle Tdec can be comparable or shorter than the already mentioned time scale w~ 1 what demands a different 
theoretical treatement [50]. 

2) Thermalization stage, when the density matrix of the Brownian particle written in momentum representation is 
close to diagonal i.e. Pt(p,p') differs essentially from zero for |p — p'| < \]1MT. The detailed analysis of this regime 
and the manifest standard form of MME which depends on the so-called dynamic structure factor characterizing the 
environment are presented in [51]. 

Another completely different situation where the continuous spectrum of Hs appears is the theory of many-body 
open systems. An ensemble of interacting quantum spins coupled to an infinite heat bath and described in the 
thermodynamic limit is a perfect example. The mathematical formalism used in this context to construct proper 
MME, difficulties and partial solutions of the problems are discussed in [52]. 
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